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THERMAL  EXPANSION  OF  ELASTIC-PLASTIC 
COMPOSITE  MATERIALS 
by 


George  J.  Dvorak 
Department  of  Civil  Engineering 
Rensselaer  Polytechnic  Institute 
Troy,  NY  12180 

Abstract 

Exact  relationships  are  derived  between  instantaneous  overall 
thermal  stress  or  strain  vectors  and  instantaneous  overall  mechanical 
stiffness  or  compliance,  for  two  binary  composite  systems  in  which  one 
of  the  phases  may  deform  plastically.  Also,  the  local  instantaneous 
thermal  strain  and  stress  concentration  factors  are  related  in  an  exact 
way  to  the  corresponding  mechanical  concentration  factors.  The  results 
depend  on  instantaneous  thermoelastic  constants  and  volume  fractions  of 
the  phases.  They  are  found  for  fibrous  composites  with  two  distinct 
elastically  isotropic  or  transversely  isotropic  phases,  and  for  any 
binary  composite  with  elastically  isotropic  phases.  The  results  indicate 
that  in  the  plastic  range  the  thermal  and  mechanical  loading  effects  are 
coupled  even  if  the  phase  properties  do  not  depend  on  changes  in  temper¬ 
ature.  The  derivation  is  based  on  a  novel  decomposition  procedure  which 
indicates  that  spatially  uniform  elastic  strain  fields  can  be  created  In 
certain  heterogeneous  media  by  superposition  of  uniform  phase  eigen- 
strains  with  local  strains  caused  by  piecewise  uniform  stress  fields 
which  are  In  equilibrium  with  prescribed  surface  tractions.  The  method 
^s  extended  U  discretized  microstructures,  and  also  to  analysis  of 
moisture  absorption  and  phase  transformation  effects  on  overall  response 
and  on  local  fields  In  the  two  composite  materials.  ^ 


1.  INTRODUCTION 


The  response  of  elastic  composite  materials  to  spatially  uniform 
changes  in  temperature  is  well  understood.  An  essential  contribution  to 
the  solution  of  this  problem  was  made  by  LEVIN  (1967),  who  found  that 
macroscopic  thermal  expansion  coefficients  of  a  composite  medium, 
consisting  of  two  distinct  isotropic  phases  of  arbitrary  shape,  depend 
In  a  unique  way  on  overall  elastic  moduli  of  the  aggregate  and  on 
thermoelastic  constants  of  the  phases.  Thus,  if  the  elastic  moduli  are 
known,  the  thermal  expansion  coefficients  can  be  calculated.  This  line 
of  inquiry  was  extended  by  SHAPERY  (1968),  who  derived  bounds  on  thermal 
expansion  coefficients  of  multi -phase  composites  with  isotropic  phases, 
while  ROSEN  and  HASHIN  (1970)  applied  LEVIN'S  approach  to  binary 
composites  consisting  of  anisotropic  phases,  and  they  also  found  bounds 
on  overall  thermal  expansion  coefficients  of  multiphase  materials. 
BUOIANSKY  (1970)  gave  self-consistent  estimates  of  several  thermal  and 
thermoelastic  properties  of  multiphase  isotropic  mixtures.  Among  the 
more  recent  contributions  to  the  subject  are  the  papers  by  LAWS  (1973) 
and  CRAFT  CHRISTENSEN  (1980). 

The  response  of  elastic-plastic  composite  materials  to  uniform 
thermal  changes  has  been  explored  only  to  a  limited  extent.  This  is  a 
more  difficult  problem  because  at  least  one  of  the  phases  is  inelastic, 
and  the  deformation  of  the  phases  and  of  the  composite  is  affected  both 
by  the  overall  thermal  change  and  by  the  current  macroscopic  mechanical 
loading.  In  the  plastic  state  the  thermal  and  mechanical  loading 

i  effects  are  coupled,  even  if  the  mechanical  properties  of  the  phases  do 

not  depend  on  temperature.  The  problem  is,  of  course,  nonlinear  and 

I 


must  be  solved  in  an  incremental  way.  Therefore,  the  connection  between 
the  two  effects  must  be  examined  within  a  combined  load  increment.  A 
representative  volume  element  of  the  composite  is  first  subjected  to  a 
certain  uniform  initial  macroscopic  stress  or  strain  and  to  a  uniform 
thermal  change;  then,  both  the  stress  or  strain  and  the  temperature 
experience  a  small  simultaneous  change  to  another  uniform  state.  Over¬ 
all  instantaneous  stiffness  and  compliance,  and  thermal  stress  and 
strain  vectors  are  sought. 

Earlier  solutions  of  problems  of  this  kind  have  been  limited  to 
simple  loading  situations  in  fibrous  composites,  such  as  pure  thermal 
change  (OE  SILVA  and  CHADWICK  1969)  or  thermal  change  combined  with 
axi symmetric  mechanical  loads  (DVORAK  and  RAO  1976).  More  recently 
(DVORAK  1983)  it  was  shown  that  the  total  overall  strain  increment 
caused  in  a  prestressed  fibrous  composite  by  a  small  uniform  thermal 
change  can  be  related  in  an  exact  way  to  thermoelastic  constants  of  the 
phases  and  to  instantaneous  overall  compliance.  No  restrictions  need  to 
be  imposed  on  the  type  of  prestress  or  on  the  matrix  constitutive  law 
except  for  plastic  incompressibility,  but  the  fiber  must  be  isotropic  or 
transversely  isotropic  and  remain  elastic.  This  result  has  been  applied 
in  analysis  of  a  composite  cylinder  element  (DVORAK  and  WUNG  1984) 
subjected  to  axi symmetric  mechanical  loading,  uniform  thermal  changes, 
and  variations  in  matrix  yield  stress. 

The  present  paper  develops  the  connections  between  overall  instan¬ 
taneous  mechanical  and  thermal  properties  in  a  more  general  way.  First, 
it  is  shown  that  the  overall  thermal  stress  and  strain  vectors  for  an 
elastic  fibrous  composite  with  transversely  isotropic  phases  can  be 
obtained  through  superposition  of  certain  uniform  fields  in  the  phases. 


and  local  fields  caused  by  a  uniform  overall  stress  or  strain.  A 
similar  result  is  derived  for  any  composite  consisting  of  two  isotropic 
phases.  These  results  are  then  utilized  to  find  instantaneous  thermo¬ 
plastic  properties  and  local  fields  of  these  composite  systems  for 
simultaneous  mechanical  and  thermal  load  increments.  Extensions  of  the 
results  to  discretized  microstructures,  and  to  additional  load  effects, 
such  as  phase  transformations  and  moisture  absorption  are  discussed  as 
wel  1 . 
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2.  ELASTIC  FIBROUS  COMPOSITE 


A  binary  composite  material  consists  of  a  matrix  reinforced  by 
aligned  and  bonded  cylindrical  fibers.  Both  phases  are  assumed  to  be 
homogeneous  and  transversely  isotropic  about  the  fiber  direction  X3.  In 
the  transverse  Xjx2-plane,  the  cross  sections  and  distribution  of  the 
phases  can  be  arbitrary  providing  that  the  composite  is  statistically 
homogeneous,  transversely  isotropic,  and  free  of  voids. 

A  representative  volume  element  V  of  the  composite  is  selected  and 
subjected  to  a  certain  loading  or  deformation  history  which  is  imposed 
through  application  of  uniform  overall  stresses  o®  or  strains  eO  to  the 
surface  S  of  volume  V,  Also,  a  certain  uniform  thermal  change  has  been 
applied  such  that  the  current  temperature  in  V  is  constant  and  equal  to 
Og.  At  this  particular  point  of  the  loading  sequence  simultaneous 
increments  of  d7  and  d9,  or  dF  and  d0,  are  applied  to  V. 

%  V 

The  response  of  the  composite  to  these  load  increments  is  described 
by  constitutive  equations: 

dr  =  MdoT  mde  ,  do  »  Ldr  -  t  d9  ,  ( 1 ) 

where  M,L  are  (6x6)  overall  stiffness  and  compliance  matrices,  and  m,Ji 
are  (6x1)  overall  thermal  strain  and  stress  vectors*. 

While  M  and  L  are  known,  we  wish  to  determine  the  vectors  m  and  t. 
To  this  end  it  is  necessary  to  specify  the  constitutive  equations  for 
field  averages  of  the  phases: 


*We  use  the  customary  notation  (HILL  1963,  LAWS  1973)  where,  except  as 
noted,  (6x6)  matrices  are  denoted  by  lightface  uppercase  Latin  letters, 
and  (6x1)  vectors  by  boldface  lowercase  Latin  or  Greek  letters.  Top 
bars  denote  overall  volume  averages. 


=  Mp  dap  +  (Hr  d9 ,  dap  =  Lp  dtp  -  ?.pd9  (r  =  f,m)  (2) 


which  are  analogous  to  (1);  f.m  indicate  the  "fiber"  and  "matrix" 
phases.  In  elastic  composites,  these  phases  are  interchangeable  and  f,m 
are  used  merely  for  convenience  of  notation. 

Since  both  the  composite  and  each  of  the  phases  are  transversely 
isotropic  about  X3,  it  is  possible  to  write  a  subset  of  (1)  and  (2) 
which  relates  the  first  two  stress  and  strain  invariants.  With  top  bars 
and  subscripts  r,f,m  omitted  in  (1)  and  (2)  one  obtains  (DVORAK  and 
8AHEI-EL-0IN  1979): 
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where  k,4,n  are  Hill's  (1964)  elastic  moduli,  E  =  n  -  a  =  Eay, 

B  “  OL,  and  oj  and  ol  are  linear  coefficients  of  thermal  expansion  in 
the  transverse  plane  and  longitudinal  direction,  respectively.  The 
strain  and  stress  invariants  are  defined  as: 


dc  a  dc  +  dc 
1  11  22 

doj  *  da22) 


de 
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a  dc 

(5) 

33 

■  <1033 

(6) 

With  appropriate  values  of  elastic  moduli  and  coefficients  a,B, 
equations  (3)  to  (6)  can  be  applied  either  to  the  composite  medium  or  to 


each  of  the  two  phases. 

As  long  as  M,L  do  not  depend  on  de ,  m  on  da,  and  i  on  dF,  the 

thermal  and  mechanical  contributions  to  dT  and  do*  in  (1)  can  be  found 

separately  and  superimposed.  By  assumption  M  and  L  are  known,  hence  the 

first  terms  on  right-hand  side  in  (1)  are  evaluated  without  difficulty 

for  any  given  da  or  dt.  To  find  m  and.t,  and  the  second  terms  in  (1), 

we  utilize  the  decomposition  procedure  of  DVORAK  (1983). 

In  the  first  step  of  the  procedure  the  fiber  and  matrix  phases  are 

separated  and  surface  tractions  which  preserve  the  current  local 
0  0 

stresses  a^  and  strains  are  applied  to  each  phase  r  >  f,m.  Alterna- 

0 

tively,  surface  displacements  corresponding  to  cp  may  be  prescribed  to 
0 

preserve  Op.  In  addition,  a  uniform  thermal  change  de  is  applied  to 
both  phases.  The  local  strains  caused  by  de  would  make  the  phases  in¬ 
compatible  if  the  composite  was  to  be  reassembled.  Therefore,  uniform 

A  A 

tractions  doj,  da^  of  as  yet  unknown  magnitude  are  applied  to  th''  phases 
simultaneously  with  de.  (The  top  hats  indicate  auxilliary  uniform 
fields  used  in  the  decomposition  and  reassembly  of  the  composite.)  This 
leads  to  the  following  uniform  strain  increments  in  the  separated 


phases: 
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In  the  second  step  of  the  procedure,  the  tractions  doj’  and  do^  must 
be  adjusted  to  assure  compatibility  of  the  phases  and  equilibrium  of 
these  tractions  at  phase  interfaces  and  on  the  surface  S  of  the  represen 
tative  volume  V.  Compatibility  and  equilibrium  require  that 


'm  f 

de  =  de 
1  1 


"m  '  ^ 
de  =  de 
2  2 


*  f  "m 

da  =  da  =  dSr 
1  I  ' 


I  111 

C  da  +  C  da  =  dSa 
f  2  m  2  ^ 


where  dSy  and  dS;\  are  surface  tractions  which  need  to  be  added  at  S  to 

preserve  overall  equilibrium  of  V  while  do'*  and  da'*  are  applied  to  the 

1  2 

phases.  The  magnitudes  of  phase  volume  fractions  c  +  c  *  1  need  also 

f  m 

be  known  at  this  point. 

All  strain  and  stress  Increments  in  (7)  are  uniform,  hence  equa¬ 
tions  (8)  to  (10)  are  exact  for  any  transverse  plane  geometry.  These 
relations  suggest  that  spatially  uniform  strain  fields  can  be  created  in 
certain  heterogeneous  media  by  superposition  of  uniform  eigenstrains 
afde,  fci^de  in  the  phases,  with  local  strains  caused  by  piecewise  uniform 
stress  fields  which  are  in  equilibrium  with  surface  tractions  dS/v,  dSy. 

Internal  equilibrium  and  compatibility  of  the  phases  in  V  depend 

•  r  *  r  “  r  ~  r 

only  on  the  eight  unknown  strains  and  stresses  de  ,  de  ,  do  ,  da  ,  and 

1  2  1  2 

not  on  dSAt  dSj.  Therefore,  (7),  (8)  and  (9i)  represent  a  system  of 
seven  equations  for  the  eight  unknwons.  To  find  a  solution,  it  is 
necessary  to  introduce  an  additional  constraint. 
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A  particular  choice  which  will  be  useful  in  the  sequel  is: 


■m  *  m 

do^  *  P  dOj^  , 


(11) 


where  p  ^  0  is  a  constant. 

Now,  dSj  and  dS^  follow  from  (7),  after  a  substitution  of  (9)  and 
( 10) ,  and  from  (8) : 


a,  dS_  +  a.  dS,  +  a,  de  =0 
1  T  2  A  3 


b,  dS  +  b„  dS  +  b-  de  =  0 
1  T  2  A  3 


where: 


.  _  ^f  r  Cfn  tf 


1  ^f 

^2  ’  ■  Cf  IfFf  •  ^3  -  af  -  am 


W  IE  “ 


(12) 

(13) 


1 

^f 


b3  =  -  8f  +  Bp 


and: 


dSy  a  ST  d9  ,  dSA  ■  SA  d0 


ST  »  (a2b3  -  a3b2)/(aib2  -  a2bi) 


SA  *  (a3bi  -  aib3)/(aib2  -  a2bi) 


(14) 

(15) 

(16) 


At  this  point  we  change  from  the  invariants  (5),  (6)  to  the  (6x1) 
stress  and  strain  vectors  and  write 


I 


9 


10 


In  the  final  step  of  the  decomposition  procedure,  the  composite  is 
reassentled  and  the  surface  tractions  dSt,  dSA  removed.  Of  course,  the 
local  strains  and  stresses  (17)  to  (23)  already  assure  that  the  phases 
are  compatible  and  in  internal  equilibrium,  in  fact  they  are  equal  to 
local  fields  caused  in  the  composite  by  simultaneous  application  of  de , 
dSA  and  dSj.  They  must  be  now  added  to  local  fields  caused  in  the 
composite  by  surface  tractions  -dS^.  -dSA« 

The  final  results  assume  a  concise  form  with  the  definitions 

h  =  [hi  hj  h2  0  0  0]^ 


Sa  “  [sj  ST  SA  0  0  O]"*" 

(24) 

Y  =  [1  1  Y  0  0  0]T 
% 

p  =  [1  1  p  0  0  O]’’’ 

where  [  denotes  a  transpose  and  the 

coefficients  appear  in  (11), 

(15),  (16),  and  (19)  to  (21). 

Therefore,  for  do  »  0,  de  *  0  in 

dl): 

de  =  m  de 

(25) 

ffl  »  h  -  MSj 

(26) 

where  m  is  the  overall  thermal  strain  vector  and  M  is  the  known  overall 
compl 1 ance. 

Also,  suppose  that  the  local  stresses  in  the  phases  are  written  in 
terms  of  concentration  factors: 

dOp  »  Bf  do  bp  d0  ,  (r  ■  f,m)  (27) 


and  that  Bf  are  known. 

With  regard  to  (22)  and  (23)  one  obtains: 


bf  =  ST  Y  -  Bf  Sa 

■  sj  p  -  Bm  Sa 

%  ^  ^ 


(28) 


Similar  results  can  be  found  for  a  fully  constrained  composite 
subjected  to  a  uniform  thermal  change.  Recall  that  the  strains  (17)  and 
(18)  are  actually  equal  to  overall  strains  under  de ,  dS/^,  dSy.  This 
follows  immediately  from  (8)  and  (21).  These  overall  strains  must  be 
removed,  and  the  local  fields  adjusted  accordingly. 

Therefore,  for  dF  =  0,  de  ^  0  in  (I2): 

do  =  -  A  de 

(29) 

A  =  -Sa  +  L  h 


where  A  is  the  overall  thermal  stress  vector,  and  L  Is  the  known  overall 
% 

sti ffness. 

Also,  if  the  local  strains  are  written  in  terms  of  concentration 
factors: 


dcp  a  Ap  dF  -  ap  de  ,  (r  =  f,m)  (30) 

and  if  An,,  Af  are  known,  one  obtains  with  the  help  of  (17)  to  (21): 

af  =■  (Af  -  I)  h  ,  an,  a  (Am  -  I)  h  (31) 

^  ^  ^  ^ 

To  facilitate  applications  we  note  that 


I*  ■*  •’  • 
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in  =  [aj  aj  (XL  0  0  0]^  (32) 

where  ay,  ol  are  linear  coefficients  of  thermal  expansion  in  the  trans¬ 
verse  plane,  and  in  longitudinal  direction. 

For  any  binary  fibrous  composite  with  known  phase  preperties  and 
phase  volume  fractions,  the  effect  of  thermal  change  is  reduced  to 
equivalent  mechanical  loads  and  to  certain  uniform  fields  in  the  phases. 
Thus  m  and  l  are  found  in  terms  of  M  and  L,  and  a-,  bp  in  terms  of  Ap, 
Bp.  All  these  relations  are  exact.  While  the  constant  p  is  a  free 
parameter,  none  of  the  results  actually  depend  on  p.  For  each  p  one 
obtains  by  superposition  a  solution  to  the  same  boundary  value  problem. 
According  to  the  uniqueness  theorem  in  the  theory  of  elasticity,  all 
such  solutions  must  coincide.  This  can  be  verified  by  numerical 
calculations. 


13 


3.  TWO  ISOTROPIC  ELASTIC  PHASES 


Suppose  that  a  composite  aggregate  consists  of  two  perfectly  bonded 
elastic  phases,  which  are  distinct  but  isotropic.  The  microstructural 
geometry  can  be  arbitrary,  providing  that  the  composite  is  statistically 
homogeneous  and  free  of  voids.  The  composite  itself  need  not  be 
isotropic,  it  can  be  reinforced  by  aligned,  braided,  or  otherwise 
distributed  continuous  fibers,  short  fibers,  particles  of  any  shape, 
and  by  combinations  of  such  reinforcements. 

Assume  that  the  overall  constitutive  relations  are  again  given  by 
(1)  and  that  the  overall  compliance  M  and  stiffness  L  are  known.  Also, 
let  local  fields  be  described  by  (27)  and  (30),  and  assume  that  Ap,  Bp 
are  known.  As  in  the  previous  section  we  utilize  the  subscripts  f,m  to 
identify  the  two  phases,  even  though  we  no  longer  require  that  either 
phase  be  of  cylindrical  shape.  Local  phase  properties  need  be  known 
only  in  terms  of  bulk  moduli  Kf,  K^,  and  linear  thermal  expansion  co¬ 
efficients  Of, 

We  again  pose  the  problem  described  in  the  previous  section:  The 

composite  has  been  loaded  by  a  certain  uniform  overall  stress  a  ,  or 

_o  _ 

strain  e  ,  and  uniform  temperature  Og.  Simultaneous  increments  do  and 

de ,  or  dT  and  de  are  applied.  The  response  of  the  aggregate  is  sought 
% 

in  terms  of  the  overall  thermal  strain  and  stress  vectors  m  and  t,  and 

phase  concentration  factors  ap,  bp,  (r  »  f,m). 

%  ^ 

As  in  Section  2,  the  phases  are  first  separated,  and  loaded  by  d9 

A 

and  by  certain  unknown  tractions  dop.  The  nonvanishing  stress  and 
strain  Increments  are 
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(33) 


.  f  .f 

dc  *  de  =  de  =  dSf/(3Kf)  +  af  d6 
II  22  33  ^  ^  ^ 


*  (n  'ni  , 


To  assure  equilibrium  and  compatibility: 
dSf  3  dSm  *  dS, 


f  m 
deij  *  <leij. 


and 


dS  =  S  d9  ,  s  »  -  3(of  -  an,)/(l/Kf  -  l/K^) 


(34) 


(35) 


The  composite  is  now  loaded  by  three  equal  overall  normal  stresses 
dS,  and  by  d6.  Local  strain  and  stress  fields  follow  from  (33)  and 
(34). 

Finally,  the  composite  is  reassembled  and  surface  tractions  dS  are 
removed. 

Let 

q  «  s/(3Kf)  +  Of  »  s/(3Kni)  + 


q  -  q  [1  1  1  0  0  0]T 
% 

s  >  s  [1  1  1  0  0  0]T 
% 

In  analogy  with  (25)  to  (30);  one  obtains 


(36) 
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dor  =  br  d9  (r  =  f,m) 

bf  =  (I  -  Bf)  s  ,  b|n  =  (I  -  B|f))  s 

For  dT  -  0,  de  *  0  in  (1): 

d7  =  -  a  d9 
^  % 

4  =  -  s  +  Lq 


(38) 


(39) 


■s  ^  N 


dcr  =  ar  <1®  (r  =  f,m) 

^  N 

af  “  (Af  -  I)  q  ,  ajjj  =  (Ajn  -  I)  q 


(40) 


If  the  conposite  is  nacroscopically  anisotropic,  then  M  and  L  can 
depend  on  up  to  21  elastic  constants  and 


m  =  [al  a2  a3  04  05 


(41) 


where  a  1  to  og  are  overall  linear  thermal  expansion  coefficients,  de¬ 
fined  by  (l^)  at  do  »  0.  For  a  fibrous  coitposite  which  Is  transversely 
Isotropic  one  recovers  (32).  For  an  Isotropic  coirposite  all  a  In  (32) 
become  equal  to 
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K. 


(af  -  am^,  1  1 


a  =  am 


Ik  •  Km) 


(42) 


where  K  is  the  overall  bulk  modulus. 

This  last  equation  was  derived  in  a  different  way  by  LEVIN  (1967) 
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4.  ELASTIC-PLASTIC  COMPOSITES 
(i)  Fibrous  composites 

Consider  again  the  fibrous  composite  system  of  Section  2.  Suppose 
that  the  matrix  phase  is  elastic  within  a  certain  stress  region,  but 
becomes  elastic-plastic  when  a  given  yield  condition  has  been  satisfied. 
The  fiber  remains  elastic  until  failure.  This  suggests  a  metal  matrix, 
which  is  usually  elastically  isotropic.  Thus  the  matrix  elastic  moduli 
in  (3)  and  (4)  become  related  as  follows: 


^m  ■  l^m  ■ 


"m 


=  k- 


^m  ■  '’m  "  "m/'^m 


Jim-  , 


(43) 


am  ® 


where  Em*  '^m  t^’®  isotropic  constants,  and  Bm  i*  the  linear  thermal 

expansion  coefficient  of  the  matrix. 

In  the  plastic  range  the  matrix  response  is  assumed  to  be  the 
piecewise  linear  and  given  by  (2),  but  M^  and  Lq,,  %  and  tpi  4^*®  f’ow 

N  > 

instantaneous  compliance  and  stiffness  matrices,  and  thermal  vectors,  at 
a  particular  point  of  a  loading  path.  We  assume  that  Mm  and  Lm  are 
symmetric,  satisfy  the  requirement  of  plastic  incompressibility  of  the 
matrix,  and  do  not  depend  on  d9;  %  and  tm  not  depend  on  de,  but  may 
be  functions  of  Bq* 

Furthermore,  we  assume  that  the  response  of  the  composite  to  any 

purely  mechanical  loading  by  uniform  do  or  by  dT  Is  also  piecewise  linear 

%  % 

and  described  by  (1).  Suppose  that  instantaneous  overall  properties  M 
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and  L  in  (1),  as  well  as  instantaneous  phase  concentration  factors  ^r, 
in  (27)  can  be  evaluated  for  any  given  mechanical  loading  step.  The  in¬ 
stantaneous  thermal  properties  m,  i,  and  concentration  factors  ar  and  bp 
are  to  be  determined. 

Therefore,  we  again  pose  the  problem  stated  in  Section  2:  A  repre¬ 
sentative  volume  V  of  the  composite  has  been  subjected  to  a  certain 

loading  or  deformation  history  such  that  the  current  overall  stresses  and 

-0  _D 

strains  in  V  are  uniform  and  have  magnitudes  a  and  e  .  Also,  the  cur- 

rent  temperature  in  V  is  constant  and  equal  to  Sq.  At  this  particular 

point  of  the  loading  sequence  we  apply  simultaneous  increments  of  da  and 

db ,  or  dT  and  de  in  V,  and  wish  to  evaluate  instantaneous  values  of  M,  L, 
> 

m,  i,  and  of  the  concentration  factors  Ap,  Bp,  ap,  bp  during  the  loading 
step. 

This  problem  is  solved  by  the  decomposition  procedure  of  Section  2. 
Initially,  the  composite  is  subjected  to  the  prescribed  thermal  change 
ds  and  to  simultaneously  applied  surface  tractions  dSj  and  dS^  given  by 
(14)  to  (16).  These  thermal  and  mechanical  loads  create  local  strain 

A  ^ 

increments  dc^,  in  (18)  and  stress  increments  da^  in  (23).  In  general, 
these  strain  increments  may  be  inelastic.  However,  since  the  matrix  is 
plastically  incompressible,  it  Is  possible  to  assure  that  these  incre¬ 
ments  correspond  to  purely  elastic  deformation  in  both  phases  under  de, 
dS/y  and  dSj.  This  is  obviously  the  case  when  one  chooses  p  =  1  in  (11), 
so  that  the  stress  and  strain  increments  in  the  matrix  are  Isotropic. 

The  tractions  dS/y  and  dSy  must  now  be  removed.  This  may  lead  to 
plastic  straining  in  the  matrix,  which  corresponds  to  or  is  caused  by 
de .  Also,  if  an  overall  stress  increment  do  is  applied  simultaneously 


••  li  -.v  Vi.  V'l  V  L 


with  de  ,  then  dS^,  dSjf  and  da  must  be  added  and  applied  together. 

The  final  results  for  a  plastically  deforming  composite  can  now  be 
written  on  the  basis  of  (24)  to  (31).  With  p  =  1  in  (11)  to  (23),  we 
retain  the  definitions  (24)  of  vectors  h,  Sg  and  y  and  of  their 
components  but  replace  that  of  p  with 

1  =  [1  1  1  0  0  0]T.  (44) 

For  loading  by  do  and  de,  the  overall  strain  increment  dF  is: 

^  V 

de  =  h  d0  +  M(do  -  S3  dft).  (45) 

A  comparison  with  (1^)  again  yields  the  form  (26) 

m  »  h  -  Msa  (46) 

We  note  that  for  p  »  1  and  an  isotropic  matrix  one  obtains  from 
(19)  and  (43)  the  following  expression  for  h  in  (24): 

h  =  h  1  (47) 

^  % 

where 

h  =■  hi  =  h2  =  ST/(3Kn,)  +  Sm, 
and  Sj  is  given  by  (15). 

Therefore,  the  first  term  in  (45)  is  an  overall  isotropic  strain 
increment,  and,  according  to  (18),  it  is  equal  to  the  matrix  strain  in¬ 
crement.  The  loading  vector  da  -  s^  ds  represents  total  mechanical  load 
that  must  be  applied  to  the  composite  to  reflect  the  effect  of  simultan¬ 
eous  application  of  da  and  d0,  M  is  the  instantaneous  overall 
compliance  corresponding  to  this  loading  vector. 

The  stresses  caused  in  the  phases  by  simultaneous  application  of  da 


and  de  are; 


duf  =  sj  Y  d0  +■  Bf(d(J  -  Sj  d6 ) 


do^  =  ST  1  d0  +  Bm(do  -  Sa  d0) 


where  Bf,  3^  are  the  instantaneous  stress  concentration  factors  for  the 
overall  mechanical  load  increment  d7  -  Sj  de.  If  these  increments  are 
described  by  (27),  then  the  instantaneous  thermal  stress  concentration 
factors  become: 


bf  ST  Y  -  Bf  Si 


bm  =  ST  1  -  Bm  Sa 


Equations  (45)  to  (49)  convert  in  an  exact  way  the  thermomechanical 
problem  into  a  mechanical  loading  problem  along  the  incremental  path 
d7  -  Sa  d9 . 

Next,  consider  loading  by  and  d9.  The  composite  is  first  sub- 
jected  to  loading  by  d0,  dS;^,  and  dST.  which  causes  isotropic  strains 
h  d0  in  both  the  composite  and  matrix.  Since  the  overall  strain  incre- 
ment  Is  now  prescribed,  the  h  de  and  any  additional  overall  strains  must 
be  equal  to  dT.  Hence,  the  overall  stress  increment  Is: 


da  a  Sa  d0  +  L(de  -  h  d9) 


A  comparison  with  (I2)  again  yields  (29) 
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While  Sj,  is  not  isotropic,  together  with  de  it  causes  an  isotropic  stress 
increment  sj  1  dO  in  the  matrix.  This  is  found  from  (23)  at  o  =  1.  Ac¬ 
cordingly,  plastic  loading  of  the  composite  is  caused  only  by  the  second 
term  in  (50).  The  overall  mechanical  strain  is  equal  to  dT  -  h  d9 .  L 
is  the  instantaneous  overall  stiffness  corresponding  to  this  strain 
i ncrement. 

The  strain  increments  in  the  phases  are: 


dcf  =  h  d9  +  Af(de  -  hd9) 


de^  =  h  d9  +  An,(de  -  hd6) 


(52) 


where  Af,  A^  are  instantaneous  strain  concentration  factors  for  overall 

strain  (k  -  hde.  If  (30)  is  used,  then  the  instantaneous  thermal  strain 
%  % 

concentration  factors  are: 


af  =  (Af  -  I)  h 
%  % 


^m  ®  (^m  ”  1) 


(53) 


Inasmuch  as  the  instantaneous  M  and  L  may  have  as  many  as  21  inde¬ 
pendent  coefficients,  the  vectors  m  in  (46)  and  i  in  (51)  may  have  6  in¬ 
dependent  coefficients.  For  example,  m  assumes  the  form  (41),  with 
(i  =  1  to  6)  representing  instantaneous  thermal  expansion  coefficients 
of  the  composite. 

As  in  the  previous  case,  (50)  to  (52)  convert  in  an  exact  way  the 


thermomechanical  problem  into  a  mechanical  deformation  problem  along  the 
path  de  -  hd6.  It  is  seen  that  in  both  cases  the  thermal  and  mechanical 
effects  are  coupled,  even  though  phase  mechanical  properties  do  not 


depend  directly  on  temperature. 


( ii )  Two-phase  composites 

Finally,  we  consider  the  two-phase  composite  with  isotropic  phases 
of  arbitrary  geometry.  Section  3.  The  reinforcement  phase  (f)  is  as¬ 
sumed  to  remain  elastic,  while  the  matrix  phase  (m)  may  become  elastic- 
plastic  when  a  given  yield  condition  has  been  satisfied.  In  the  plastic 
region,  the  ratrix  constitutive  relation  is  described  by  (2),  with  Mn, 
and  Lm  replaced  by  instantaneous  compliance  and  stiffness.  Again, 
and  are  assumed  to  be  piecewise  linear,  symmetric,  and  satisfy  the 
requirement  of  plastic  incompressibility  of  the  matrix.  Also,  we  assume 
that  overall  instantaneous  properties  M  and  L  of  the  composite,  as  well 
as  the  instantaneous  concentration  factors  A^,  can  be  evaluated  for 
any  purely  mechanical  overall  stress  or  strain  increment  in  the  elastic 
and  plastic  range. 

To  find  instantaneous  thermal  properties  m,  t,  and  the  concentre- 

tion  factors  a^,  b^,  we  again  consider  a  representative  volume  V  of  the 

_o 

composite  which  has  been  loaded  to  current  uniform  overall  stress  a  . 

—0 

strain  z  ,  and  temperature  Bq*  The  volume  V  is  now  subjected  to  addi¬ 
tional  increments  of  d7  and  de,  or  dT  and  de .  As  in  Section  3,  we  apply 
overall  increments  of  temperature  de  and  of  isotropic  stress  dS,  with  dS 
given  by  (35).  Resulting  phase  stresses  and  strains  follow  from  (33), 
they  are  isotropic  and  by  assumption  cannot  cause  plastic  deformation  in 
the  matrix.  The  surface  tractions  or  strains  must  now  be  adjusted  to 
satisfy  the  prescribed  d7  or  dF  at  the  boundary  S  of  V. 

For  the  case  of  de  and  d7  applied  simultaneously  one  obtains  the 
overall  composite  strain  Increment 


de  =  q  d0  -t-  M(do  -  s  d0) 


(54) 


where  q  and  s  are  given  by  (35)  and  (36),  and  M  i s  the  overall  instan¬ 
taneous  compliance  for  the  mechanical  stress  increment  do  -  s  de. 

A  comparison  with  (1^)  again  yields  the  form  (37) 

m  =  q  -  Ms  (55) 

The  stresses  caused  in  the  phases  are: 

dof  =  s  d0  +  8f(d7  -  s  de) 

(56) 

don,  =  s  d0  t-  Bm(do  -  s  de) 

where  Bf,  3^  are  instantaneous  stress  concentration  factors  for  the 
overall  mechanical  load  increment  do  -  s  d0.  From  (27)  and  (56),  the 
instantaneous  thermal  stress  concentration  factors  are: 

bf  =  (I  -  8f)s  ,  bm  =  (I  >  Bn,)s  (57) 

For  the  case  of  d9  and  dc  applied  together,  one  obtains  the  overall 
stress  increment 


do  =  s  de  +•  L(d€’  -  q  do) 


(58) 


and 


t  =  -  s  +  Lq 


(59) 


where  L  is  the  instantaneous  composite  stiffness  for  the  overall  strain 
increment  ck  -  q  d9 ,  and  q  is  given  by  (36).  The  phase  strain  increments 
in  this  case  are 
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de|r  =  q  d9  +  Af(de  -  q  d0 ) 


=  q  d9  +  An,(dr  -  q  d0) 

and  the  instantaneous  concentration  factors: 


(60) 


af  =  (Af  -  I)q  ,  a,ii  =  (A,^  -  I)q  (61) 

^  ^  X 

where  Af,  A^  are  instantaneous  strain  concentration  factors  for  an  over 
all  mechanical  strain  increment  equal  to  dT  -  q  d9. 

Equations  (54),  or  (58),  again  convert  in  an  exact  way  the  thermo- 
mechanical  loading  problem  into  a  mechanical  one  along  a  loading  path 
d7  -  s  d9 ,  or  a  strain  path  dF  -  q  d0 ,  respectively.  As  in  the  case  of 
a  fibrous  composite,  the  thermal  and  mechanical  loading  effects  are 
coupled  in  the  instantaneous  M,  L,  and  Ap,  Bp. 


5.  DISCRETIZATION  OF  THE  PHASES 


Results  of  the  previous  sections  depend  on  the  availability  of 
overall  instantaneous  mechanical  stiffnesses,  compliances,  and  phase 
concentration  factors  in  each  loading  step.  These  quantities  need  to  be 
evaluated  for  a  certain  model  geometry  of  the  composite  material.  An 
important  consideration  in  the  choice  of  a  material  model  is  the  fact 
that  the  thermal  loading  paths  do  -  s^  de  in  (45)  and  do  -  s  de  in  (54), 
as  well  as  the  thermal  strain  paths  (f  -  h  de  in  (50)  and  (C  -  q  de  in 
(58),  may  have  a  significant  isotropic  component.  That  is  easily  seen 
from  the  definitions  (24),  (36)  and  (47)  of  s^,  s,  h  and  q.  It  follows 
that  the  material  model  chosen  for  analysis  of  the  mechanical  response 
must  give  reasonably  accurate  predictions  when  the  composite  is  loaded 
by  isotropic  overall  stresses  or  strains.  This  restriction  may  exclude 
certain  models  which  are  primarily  useful  in  predicting  the  behavior  of 
a  fibrous  lamina  under  in-plane  loads,  such  as  the  VFO  model  (DVORAK  and 
BAHEI-EL-OIN  1982). 

Another  important  consideration  in  the  choice  of  a  material  model 
is  the  fact  that  when  the  matrix  phase  becomes  plastic,  the  local  pro¬ 
perties  (2)  are  stress-dependent,  and  therefore,  L^,  and  are  no  longer 
spatially  uniform.  Even  If  (2)  are  regarded  as  relations  for  averages 
In  the  phases,  the  phase  properties  need  to  be  determined  for  the  actual 
local  fields  or  their  approximations.  This  excludes  application  of 
certain  averaging  techniques,  such  as  the  self-consistent  method,  which 
assume  that  phase  fields  are  uniform. 

These  considerations  suggest  that  the  chosen  composite  model  should 
be  based  on  a  specific  representative  geometry  of  the  microstructure. 
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which  allows  for  discretization  of  each  phase  into  a  nunter  of  sub¬ 
elements  with  locally  uniform  fields.  An  example  of  such  an  approach 
was  outlined  by  DVORAK  and  TEPLY  (1985).  In  general,  if  the  represen¬ 
tative  volume  of  the  composite  and  the  phase  geometry  in  this  volume  are 
specified,  then  each  phase  can  be  subdivided  into  a  certain  nunfcer  of 
finite  elements,  and  the  overall  properties  L,  M,  local  properties  in 
the  plastically  deforming  subelements,  as  well  as  the  concentration 
factors  for  each  subelement  can  be  calculated  for  any  load  or  strain 
i ncrement. 

Suppose  that  the  subelement  stresses,  strains,  as  well  as  the 
stiffnesses  and  compliances  of  plastically  deformed  subelements  have 
been  found  for  a  certain  increment  d7  or  dT  applied  to  the  representa- 
tive  volume  at  d9  =  0.  Let  subscripts  i,  and  j,  denote  subelements  in 
the  matrix,  and  fiber,  respectively.  If  the  partial  contributions  of 
ea  ,h  row  of  d7  or  d?  are  identified,  one  can  write  the  uniform  subele¬ 
ment  fields  in  the  form 


^oitn  =  Bitn  ‘lojf  =  Bjf  do 

^‘^im  “  d7  de  jf  =  Ajf  dr 

^  ^  ^  % 


(62) 


where  the  A,  B  are  instantaneous  subelement  concentration  factors. 

One  can  also  write  the  following  relations  between  the  overall 
averages  and  the  uniform  local  fields  In  the  subelements  of  the  repesen 
tative  volume: 


do  »  I  c^  doin,  +  E  Cj  dojf 
^  ^  > 

(63) 


de  •  Z  C^  deiffl  ♦  E  Cj  dcjf 
^  ^  ^ 
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where  ,  cj  are  subelement  volume  fractions  such  that 

E  Ci  =  Cn,  ,  j:  Cj  =  Cf  ,  Cf  +  Cn,  =  1  (64) 

Using  (62)  and  (63)  one  can  obtain  the  average  phase  concentration 
factors 

^  ^1  ^im  ®m  “  cT  ^  ^1  ®1m 

(65) 

Af  =  I7  ^  <^1  Ajf 

From  the  local  Instantaneous  In  the  plastically  deforming 

subelements,  known  L^,,  f>n.  Lf,  Mf  In  the  elastic  subelements,  and  (62) 
to  (65),  one  can  find  the  overall  instantaneous  properties  as: 


2  Ci  Lim  Lf  Ajf 

M  Ci  M^m  8im  +  E  Cj  Mf  Bjf 


(66) 


The  local  thermal  strain  vectors  1%  and  mf  remain  constant  In  each  sub- 
element,  at  least  for  a  given  de,  and  equal  to  those  of  the  elastic 
phase.  These  thermal  stress  vectors  are: 


^1m  =■  Mm  mn,  ,  f.f  »  Lf  mf  (67) 

The  decomposition  procedure  can  now  be  applied  to  the  discretized 
representative  volume.  The  results  follow  from  those  presented  In 
Section  4. 

In  the  fibrous  composite  one  obtains;  For  do  ^  0,  de  ^  0,  the 
overall  strain  Increment  dF,  and  overall  m,  follow  from  (45)  and  (46), 
with  M  taken  from  (66).  The  local  subelement  stresses  and  thermal 
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stress  concentration  factors  are,  in  analogy  with  (48)  and  (49): 


dojf  =  sj  Y  d0  +  Bjf(do  -  Sj  d6 ) 


doim  =  ST  1  d9  +  8im(do  -  Sj  do ) 


bjf  =  ST  Y  -  Bjf  Sa 


(68) 


.■ 

>.v  • 


•:>r: 


*5im  =■  ST  Y  -  Qim  Sa 
^  ^  ^ 

For  dT  ^  0  do  ^0,  the  overall  stress  increment  da  and  overall  t, 
follow  from  (50),  (51)  and  (66),  and  the  local  fields  and  thermal  strain 
concentration  factors  are  as  in  (52),  (53): 

dcjT  a  h  do  +  Ajf(de’  -  h  do) 


dcin,  =  h  do  +  Ain,(de  -  h  do ) 

'  '  (69) 

^jf  “  (Ajf  -  0  ^ 


^im  *  (Aim  -  I)  h 

These  results  can  be  utilized  to  find  average  instantaneous  thermal 
strain  concentration  factors  in  the  phases  as: 


^  ^  Cl  Urn 


*  c7^ 

*=1  .^im 


(70) 
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6.  RELATED  APPLICATIONS 


In  addition  to  mechanical  and  thermal  loading,  the  composite  may 
also  undergo  a  phase  transformation  such  that  one  or  both  phases,  if 
free,  experience  a  volume  change  over  an  increment  of  temperature  de: 

dVr/V  =  deJi^  =  3 

Furthermore,  composites  with  polymer  matrices  may  absorb  moisture.  This 
causes  swelling  of  the  unconstrained  matrix  material.  If  the  moisture 
concentration  is  uniform,  then 

dekk  *  3  Yn,  dc  (72) 

where  yri  the  linear  swelling  coefficient  and  c  is  moisture  concen¬ 
tration.  If  the  matrix  remains  elastic  in  dilatation,  and  inviscid, 
then  the  above  theory  may  be  applied  with  the  following  adjustments. 

Suppose  that  de  and  dc  are  applied  simultaneously  and  that  a  free 
phase  r  undergoes  total  volume  change 

p  9  T  H 

=  ^<^kk  <^kk  *  <kkk)  (73) 

r 

If  this  superposition  holds  in  a  piecewise  linear  manner,  then  one 
can  write  in  phase  r  for  each  loading  step 


r  r  r 

dckk  =  3(ap  +  Yr)  *  (74) 

where,  in  a  particular  loading  step  n; 

r  9  T  r  0  H 

5n  »  (<kkk/<lekk)^  .  Cn  *  (<t kk/<lEkk)^ 

are  known  distribution  coefficients.  Hence. 


dtid^  =  3  (1  +  fin  +  Cn)  '^r  <1®* 


(75) 


This  suggests  that  for  each  loading  step  n  one  can  evaluate  a 
certain  multiplier  of  3de  that  can  be  substituted  for  the  Instantaneous 
linear  thermal  expansion  coefficient  of  phase  r  In  (2).  Indeed,  even  In 
the  case  of  thermal  loading  alone  It  may  be  appropriate  to  change 
with  temperature,  and  that  Is  obviously  possible  In  the  present  theory. 


7.  DISCUSSION 


Although  the  results  are  valid  only  for  the  two  binary  systems, 

they  apply  to  most  composite  materials  of  practical  Interest.  For  the 

composite  systems  In  question,  the  elastic  values  of  m  and  i  found  from 

%  % 

(26)  and  (29)  are  Identical  with  those  that  can  be  calculated  from  LEVIN 
(1967)  formulae,  or  equation  (2.20)  In  ROSEN  and  HASHIN  (1970)  and 
equation  (33)  In  LAWS  (1973).  However,  the  methods  used  In  deriving 
these  respective  equations,  and  their  Internal  structure,  are  entirely 
different.  The  decomposition  used  herein  makes  It  possible  to  find 
overall  thermomechanical  response  of  the  composite  In  the  plastic  range 
In  terms  of  Instantaneous  overall  mechanical  properties  and  thermo- 
elastic  constants  of  the  phases.  Also,  average  Instantaneous  phase 
stresses  and  strains  are  found  In  terms  of  mechanical  concentration 
factors.  All  these  relations  are  exact.  They  make  It  possible  to  con¬ 
vert  any  available  facility  for  analysis  of  Isothermal  elastic-plastic 
behavior  of  the  two  composite  systems  to  one  which  can  analyze  the 
effect  of  both  mechanical  loads  and  uniform  changes  In  temperature,  as 
well  as  other  transformation  strains  In  the  phases. 

It  has  not  escaped  our  attention  that  equations  similar  to  (7)  to 
(10)  can  be  written  for  three  phase  fibrous  composites.  The  resulting 
system  has  at  most  one  solution,  and  If  It  exists  It  leads  to  evaluation 
of  overall  thermal  vectors  and  local  thermal  concentration  factors  for 
the  elastic  three  phase  materials  which  are  analogous  to  those  derived 
In  Section  2.  However,  no  additional  constraints  are  allowed  In  this 
case  hence  (11)  cannot  be  Introduced,  and,  therefore.  It  Is  not  possible 
to  analyse  elastic-plastic  deformation  of  the  three  phase  aggregate  by  the 


method  of  Section  4(i) 
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